Zeeman levels with exotic field dependence 
in the high field phase of a S = 1 Heisenberg antiferromagnetic chain 
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We have performed electron spin resonance measurements over a wide frequency and magnetic 
field range on a single crystal of the S = 1 quasi-one-dimensional Heisenberg antiferromagnet 
Ni(C5Hi4N2)2N3(PFe). We observed gapped excitation branches above the critical field (H c ) where 
the Haldane gap closes. These branches are analyzed by a phenomenological field-theory using 
complex-field (f> 4 model. A satisfactory agreement between experiment and theory is obtained. 
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Recently, field-induced phenomena in quantum spin 
systems have attracted considerable interest. These 
include magnetization plateaus El EL Hi an d field- 
induced long-range order (LRO) 0,13,131. Gapped one- 
dimensional (ID) spin systems subject to an external 
magnetic field strong enough to close the gap are driven 
into a new phase. When the system has XY or Heisen- 
berg symmetry, this phase is critical and its low-energy 
physics is described by the Tomonaga-Luttinger liq- 
uid Q, whose elementary excitations are of the particle- 
hole type (spinon pairs). The high-energy physics which 
cannot be described using this picture has been recently 
investigated experimentally 0, 0, El as we U as theoreti- 



cally Uil. 

One-dimensional Heisenberg antiferromagnets (HAFs) 
with integer spin are typically gapped and have a disor- 
dered spin liquid singlet ground state [12|. The first ex- 
cited state is a triplet and one of the Zeeman-split triplet 
branches goes down in magnetic field. At some criti- 
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FIG. 1: The magnetic field versus temperature phase dia- 
gram of NDMAP; solid lines are guides to the eye (figure 
from Q). The ESR measurements were performed along the 
dashed line. 



cal field H c the gap closes, a qualitatively new ground 
state emerges and magnetism recovers. Above H c , a 
LRO is expected to occur in quasi- ID integer-spin HAFs 
owing to interchain interaction. Such a field-induced 
LRO in the 5=1 quasi-lD HAF was found recently 
by specific heat measurements on a single crystal of 
Ni(C 5 Hi4N2)2N 3 (PF 6 ), abbreviated NDMAP 0,11. Fig- 
ure ^ shows the magnetic field versus temperature phase 
diagram of NDMAP. The regions below and above the 
boundary denoted as ifLRo(T) correspond to the Hal- 
dane disordered phase, and to the LRO phase, respec- 
tively. The critical field H c in the present paper is defined 
as if lro extrapolated to T = 0. 

Quite recently, three distinct excitation branches were 
observed in inelastic neutron scattering (INS) experi- 
ments on NDMAP in the LRO phase above H c for H \\ a, 
and a satisfactory desription within a phenomenological 
field theory was obtained 9] . This observation is remark- 
able because one would expect two branches from the 
conventional spin- wave theory. 

Electron spin resonance (ESR) is one of the powerful 
methods to investigate magnetic excitations, especially at 
low energies with a much better energy resolution than 
INS. In Ref.0, we observed ESR signals in a single crys- 
tal of NDMAP and interpreted the results as the coexis- 
tence of ID and 3D excitations, because no satisfactory 
description was available above H c at that time. 

In the present paper we extend the ESR study of 
NDMAP to a much wider frequency range from which 
we obtain new results. A nontrivial field dependence of 
the resonance modes, very different from that expected 
from a conventional spin-wave theory, is observed. We 
compare the experimental results with model calculations 
based on the field theory used in the analysis of the re- 
cent INS experiments |9j, and obtain a good agreement 
with the theory. 

First we summarize the crystal and magnetic proper- 
ties of NDMAP. This compound crystallizes in the or- 



2 





i i 














units) 


T=1.5 K 


H // a axis 
227.8 GHz 




ESR signal (Arb. 


... 

1 1 1 


H // c sxis 
474.7 GHz 

1 





2 4 6 8 10 1 2 14 



Magnetic field (T) 

FIG. 2: An example of ESR spectra for H \\ a and H \\ c 
obtained at 1.5 K. 



thorhombic system (space group Pnmn) [14|. The Ni 2+ 
ions are bridged by azido ions forming antifcrromagnetic 
chains along the c axis. All the Ni sites in a chain are crys- 
tallographically equivalent, so that no staggered compo- 
nents of the magnetic moments are retained in the ground 
state below H c . The physics in NDMAP is, therefore, 
essentially different from that in another typical Hal- 
dane chain material NifC^HsN^NO^fClCU), abbrevi- 
ated NENP, extensively studied by ESR |ll| These stud- 
ies showed that the Haldane gap did not close at IT due 
to a staggering of the local crystal fields for Ni 2+ [l^, lljj . 
Consequently, in NENP no clear transition from the dis- 
ordered to LRO phase is observed and the ESR transi- 
tions from the ground state are always allowed. 

From the analysis of the magnetic susceptibility data, 
the following values are obtained 0: J/k B = 30.0 K, 
D/J = 0.3, 0ii = 2.10 and g± = 2.17, where J is the 
intrachain exchange interaction constant and gu and g± 
are the g values parallel and perpendicular to the chain, 
respectively. From the INS measurements done on single 
crystals of deuterated NDMAP |l8| the following param- 
eters were determined: J = 2.85 meV(= 33.1 K), J' x = 
3.5 x 10~ 4 meV (= 4.1 x 10~ 3 K), J' y = 1.8 x 10~ 3 meV 
(= 2.1 x 10" 2 K) and D = 0.70 meV (= 8.1 K), where 
J' x and J' y are the interchain exchange interactions along 
the a and b axes, respectively. 

The single crystals of NDMAP were grown in the same 
manner as that described in RefH The ESR measure- 
ments were performed using a millimeter vector network 
analyzer with extensions up to about 700 GHz from the 
AB Millimeter, France and a superconducting magnet up 
to 16 T from the Oxford Instruments, UK. We performed 
the ESR measurements in H \\ a and H \\ c geometries at 
the lowest temperature available with our spectrometer 
(1.5 K). Faraday and Voigt configurations were used for 
H || a and H \\ c, respectively. As shown in Fig. ^ the 
external magnetic field of k 9 T for ff || a brings the 
system into the LRO phase at 1.5 K. 



Figure shows a typical example of ESR spectra ob- 
served for H || a and H \\ c. The signal at about 13.5 T 
for H || a comes from the LRO phase. In Fig.[3{a), the 
resonance frequencies for H || a as a function of magnetic 
field are shown with solid circles. We also plot with solid 
squares the gap energies obtained from the recent INS 
measurements at 30 mK 9]. In the figure, the vertical 
broken line and the vertical dotted line denote the T = 
critical field H c and the boundary field 7?lro between the 
disordered and LRO phases at 1.5 K, respectively. It is 
obvious that the slope of the branches changes at H c and 
not at i?LR.o • It should be remarked that a similar change 
of slope was observed in the quasi- ID material NENC [8j 
and in the 3D weakly coupled dimer system TICUCI3 
[ToL H^ | . The observed neutron data below H c represent 
the magnon triplet split by the single ion anisotropy in 
the Haldane phase. Below H Cl only temperature-induced 
ESR transitions between the magnon branches are al- 
lowed. Only one ESR line, corresponding to the transi- 
tion between the lowest two branches, was observed near 
the critical field below H c . Above H c , three ESR lines 
were observed. Their positions coincide with the INS 
data for the gaps at q — n, confirming that they cor- 
respond to the transitions from the ground state. It is 
worthwhile to note that these lines appear at H > H c and 
not at H > .Hlro as one would naively expect, and that 
there is no anomaly in the behavior of resonance lines 
at -Hlro Indeed, the gap of the lowest magnon branch 
at q = 7r closes at H = H c at any temperature, causing 
"condensation" of q = it magnons into the ground state 
at H > H c . The amplitude of the q = tt component be- 
comes nonzero immediately above H c , making the ESR 
transitions from the ground state to the excited ones at 
q = 7r possible. The transversal staggered LRO, how- 
ever, emerges only at the higher field Hlro, because the 
phase coherence of the condensate wave function, directly 
related to the LRO, is destroyed by thermal fluctuations 
in the field range H c < H < Hlro- When the phase 
coherence sets in at H — i?LRO ; there is no peculiarity in 
the amplitude and thus no essential change in the ESR 
signal. 

Figure Bl^b) shows a frequency versus field plot for 
H || c. As in the H || a case, we observed only 
one resonance line below H c which corresponds to the 
temperature-induced transitions between the lowest two 
magnon branches. Above H c , two ESR lines were ob- 
served. For this geometry, there is no LRO at 1.5 K, so 
that the measurements were conducted completely out- 
side the LRO phase. One of the remarkable features in 
this figure is a change in slope at about 7 T, far above 
the critical field H c ~ 3.5 T, resulting from an avoided 
crossing between the two branches (see below). 

Now, we compare the experimental results with the- 
oretical calculations using a field-theoretical model 
similar to that proposed for dimerized S = ^ chains in 
Ref. l20i It is a Ginzburg-Landau-type theory written in 



3 



terms of a complex triplet field = A + iB which is 
assumed to be small, |0| <§; 1. For a dimerized S = 5 
system, the field describes the dimer wave function 
|0) = (1 - I0I 2 ) 1 / 2 !*) + Ea=x,»,>al*a), where \s) and 
\t a ) are the singlet and three triplet states of the dimer; 
it is a continuum analog of the bond boson description 
used in |l9|. The uniform and staggered magnetization 
are given by M = (0* x 0) and L = (0* + 0)(1 - |0| 2 ), 
respectively. The effective Lagrangian density in the con- 
tinuum limit has the following form: 

L = h(A-d t B-B- d t A) - (l/2)v 2 (d x A) 2 

- ^K^ + ^} + 2Jf.(4xB) (1) 

- \(A 2 ) 2 - Ax(A 2 B 2 ) - \ 2 {A- Bf. 

This model, though derived for dimerized S = h chains, 
may also be used for other gapped ID systems which are 
in the same universality class, particularly for S = 1 Hal- 
dane chains. In this case, the quantities and A ; 

should be treated as phenomenological parameters; the 
direct microscopic justification of this model for S = 1 
chains is the subject of ongoing work. The spatial deriva- 
tives of B are omitted in Q because they appear only in 
terms which are of the fourth order in A and B. By in- 
tegrating out the "slave" B-field one obtains an effective 
real-field <^ 4 -type theory similar to that of Affleck [ilj : 
actually, Affleck's theory may be viewed as a special case 
of Q with isotropic m a — fa and with simplified inter- 
action term \\ 2 = 0. Another special case m a = 
leads to the spectra which for H < H c exactly coincide 
with those obtained in the approach of Tsvelik who 
has proposed an effective theory involving three massive 
Majorana fields, and also with the perturbative formulas 
of Refs. |U If one additionally assumes the simplified 
interaction potential Ai,2 = 0, the theory becomes equiv- 
alent to that of Mitra and Halperin |l6j who postulated 
a bosonic Lagrangian to match Tsvelik's results for the 
field dependence of the gaps below H c . 

The solid lines in Figs. Et a ) and (b) represent the 
best fit of the calculated energy gaps at q = it to the ex- 
perimental data. This fit is obtained using the following 
set of model parameters (all values in meV): m a = 0.5, 
m b = 0.8, m c = 5.29, fn a = 0.353, m b = 0.342, 
m c = 0.647, A = 1, Ai = -A 2 = 0.13. We obviously 
see a change in slope of each excited branch and re- 
opening the energy gap of the lowest branch at H c . Be- 
low H c ~ 5.5 T in Fig. Eta), three excitation branches 
are in a good agreement with the neutron data. The 
ESR transitions from the ground to the excited states 
at q = 7r are forbidden because no staggered component 
exists in the ground state of this compound. In both 
figures, the ESR resonance points only appear near H Cl 
because these correspond to the transition between the 
excited states and the intensity is related to the popula- 
tion of the excited states at 1.5 K. The thin broken line 
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FIG. 3: The ESR resonance fields are plotted with solid circles 
in the frequency versus magnetic field plane for (a) H \\ a 
(the field perpendicular to the chain) and (b) H \\ c (the field 
parallel to the chain), respectively. The broken vertical line 
is the critical field H c at T = K. The solid lines are fits 
based on the field-theoretical approach. The thin broken line 
below H c denotes transition between the lowest two excited 
states. In the upper panel, the solid squares are taken from 
the neutron inelastic measurements at the wave vector q — 
7r. The dotted vertical line denotes the boundary between the 
disordered and the LRO phases at T = 1.5 K (see Fig. 0. 



in each figure indicates the transition between the lowest 
two magnon branches at q = tt. For those temperature- 
induced transitions, the agreement between experiment 
and calculation is satisfactory for H \\ a and is not good 
for H I c. Above H c , the calculated branches are close 
to the experimental data near H c , but deviate from them 
at high fields. This is probably because the theory is ap- 
plicable only for the state with a small staggered order 
and thus only for fields which are not too far above H c . 

In these calculations, we have taken into account the 
tilting of the crystal field axis of the Ni 2+ ions from the 
crystallographic c axis. The tilting angle is 15.9°. In 
Fig. Etb), a change in slope at about 7 T observed in 
the A and B branches are reproduced in the calculation. 
Because of the tilting of the crystal field axis, the exter- 
nal magnetic field applied parallel to the c axis mixes the 
wave functions of A and B branches. The origin of this 
change in slope is different from that of the changes vis- 
ible near the critical field. The observed ESR line below 
H c in Fig. Etb) deviates from the expected thin broken 
line, although a satisfactory agreement is obtained for 
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H || a; the reason for this discrepancy is not clear at 
present. For H \\ c, we did not observe ESR signals cor- 
responding to the lowest calculated mode above H c ; we 
believe this is caused by a strong damping of the lowest 
magnon mode due to the interaction with domain walls, 
because our measurement was done at a high tempera- 
ture. A similar problem arose in INS measurements 
in the H \\ a geometry done at 2 K where the lowest 
magnon mode appeared to be quasielastic, and only af- 
ter lowering the temperature to 30 mK the mode was 
successfully observed The damping should be much 
stronger in the H \\ c geometry, due to the lower energy 
and respectively higher density of domain walls. 

Remarkably, the observed field dependence of the ESR 
modes for both geometries is very different from that ex- 
pected from a conventional spin-wave theory. Indeed, in a 
conventional (classical) antiferromagnet with easy-plane 
anisotropy D and external magnetic field H applied in 
the easy plane, there are two resonance modes with the 
energies z\ = g^BH and e 2 = 2SV2DJ. For H perpen- 
dicular to the easy plane, one would classically expect to 
see one field-dependent resonance mode with the energy 
e' 2 = \J (gfiBH) 2 + e\, the energy of the other mode be- 
ing zero. We have observed and described theoretically 
completely different behaviors above H c for both direc- 
tions: in the H || a case (in-plane geometry) the energies 
of all modes increase approximately as g/is-ff, and for 
H || c we observe the (avoided) crossing of two modes 
whose energies behave roughly as g\isH and 2g^sH. 

In conclusion, we have performed ESR measurements 
on a single crystal of the S=l quasi-one-dimensional 
Heisenberg antiferromagnet NDMAP and observed ESR 
spectra corresponding to the gapped state above the crit- 
ical field where the Haldane gap closes. For H \\ a our 
ESR data above H c , taken at T — 1.5 K, agree very 
well with those from inelastic neutron scattering ob- 
tained at much lower temperature, and the observed res- 
onances are not affected by the boundary of the long- 
range-ordered phase at H ~ 9 T. For H || c we observed a 
change in slope of the ESR line resulting from an avoided 
crossing between the two modes having nontrivial slopes 
H and 2H. The observed field dependence of the magnon 
modes for both geometries H \\ a and H || c is very dif- 
ferent from that expected from a conventional spin- wave 
theory. Model calculations based on the phenomenologi- 
cal field-theoretical approach have reproduced all the fea- 
tures of the observation and a satisfactory quantitative 
agreement between theory and experiment is obtained. 
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